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Abstract.  Asymptotic  behavior  of  solutions  of  ordinary  differential  equations  in  the 
case  where  the  right  part  is  a  product  of  two  differential  expressions  with  almost  constant 
coefficients  is  investigated. 
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1.  Introduction 

The  goal  of  this  paper  is  to  obtain  asymptotic  formulas  as  x  — )•  +oo  for  the 
fundamental  system  of  solutions  of  the  following  differential  equation: 

h{h{y))  =  >^y,  X  G  i?+ :=  [0, +oo),  (1) 

where  A  is  a  complex  parameter, 

h{y)  ■■=  +  (ai  +pi(x))y(”"^^  +  (a2  +  . . .  +  (un  +Pn{x))y,  (2) 

Uy)  :=  y^^'^  +  (6i  +  +  (62  +  g^(x))y(— 2)  +  . . .  +  (6^  +  ^^{x))y,  (3) 

and  all  derivatives  are  understood  in  the  sense  of  distributions. 

Throughout  this  paper  it  is  assumed  that  ai,  02, . . . ,  On,  61,  62, . . . ,  6m  are  com¬ 
plex  numbers,  pi,p2,  ■  ■  ■  ,Pn,  qi,q2,  ■  ■  ■  ,qm  are  complex- valued  measurable  func¬ 
tions  on  i?_|_  such  that 

n 

\pi\  +  (1  +  \P2-Pi\)  \Pj\  G  LlciR+) 

j=2 
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and 

m 

kil  +  (1  +  \q2  “  9iI)  X/ ^  ^ioc{R+)- 

i=2 

More  precise  conditions  on  the  functions  pi  (i  =  1,  2, . . . ,  n)  and  qj  (j  =  1,2, ...  ,m) 
will  be  given  in  the  formulations  of  theorems  below. 

The  novelty  of  the  obtained  theorems  in  comparison  with  the  well-known 
classical  results  is  that  we  assume  that  derivatives  are  understood  in  the  sense 
of  distributions.  In  addition,  the  results  are  new  for  the  functions  pi  and  qj, 
differentiable  a  sufficient  number  of  times. 

2.  The  product  of  quasidifferential  expressions 

Let  oi,  02, . . . ,  On  be  complex  numbers  and  let  pi,p2, . .  ■  ,Pn  be  measurable 
complex-valued  functions  on  such  that 

n 

\pi\  +  (1  +  \P2  -  pi\)'^\pj\  G  L\^c{R+)^  (4) 

j=2 

where,  as  usual,  is  the  space  of  functions  Lebesgue  integrable  on  any 

segment  [a,/3]  C  R+. 

Let  F  :=  (fij)  be  an  n-dimensional  square  matrix  of  the  form 


(  0 

1 

0 

.  .  0 

0  \ 

0 

0 

1 

.  .  0 

0 

0 

0 

0 

.  .  1 

0 

/n-1.1 

fn-1,2 

fn-1,3  ■ 

fn—l,n—l 

1 

\  fnl 

fn2 

fn3 

fn—l,n 

fnn  / 

The  elements  of  the  first  (n  —  2)  rows  of  the  matrix  F  are  determined  by  the 
equalities  fij  ;=  0,  if  j  /  i  -|-  1  and  fij+i  :=  1.  Define  the  elements  fn-i,j  of  the 
penultimate  row  of  this  matrix,  assuming 

fn—l,j  ■—  Pn+l—j:  j  —  1,2,. ..71  1,  fn—l,n  •—  1; 

and  the  elements  fnj  of  the  last  row,  assuming 

fnl  ■—  {P2  ~  Pi  0,l)Pn  Cini  fnn  ■—  P2  Pi  Q-l) 

fnj  ■ —  Pn+2—j  {p2  Pi  (^l)Pn+l—j  dn+l—j^  j  —  2,  3,  ...71  1. 
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Following  the  generally  accepted  procedure  (see,  for  example,  [1,  section  I, 
p.8]),  we  define  quasiderivatives  y^p  {j  =  0, 1,  ...n)  of  a  given  function  y  by  means 
of  the  matrix  F,  assuming 

n—1 

y^p:=y‘'^\  j  =  0,  l,...,n-2,  := 

1=1 


yP  ■=  {Vf  -'^fnjVp 
1=1 

and  the  quasidifferential  expression  Tpiy),  assuming 

Tpiy)  ■=yP- 

The  domain  D{tf)  of  expression  Tpiy)  is  the  set  of  all  complex-valued  func¬ 
tions  y  for  which  the  quasiderivatives  up  to  the  order  (n  —  1)  exist  and  are 
absolutely  continuous  on  every  segment  [a,/?]  C  R+.  It  is  obvious  that  TF{y)  G 
^Zoc(^+)  any  y  G  D{tf). 

In  [2]  it  was  shown  that  for  all  y  G  D{tf)  distribution  li{y)  (see  (2))  is  a 
regular  generalized  function  and 

h{y)=TF{y).  (5) 

Further,  let  6i,  62;  •  •  • ,  be  complex  numbers,  and  yi,  (72,  •  •  • ,  fZm  be  complex¬ 
valued  measurable  functions  on  Rp  such  that 

m 

\qi\  +  (1  +  \q2  “  9il)  X/  ^  ^lociR+)-  (6) 

1=2 

Denote  by  G  an  m-dimensional  square  matrix  of  the  same  structure  as  the 
matrix  F  but  determined  by  the  numbers  bj  and  functions  qj.  The  matrix  G 
generates  quasiderivatives  y]^  {j  =  0,1,..., m)  and  quasidifferential  expression 
Tciy)-  As  before,  for  all  y  G  D{tg)  distribution  hiy)  (see  (3))  is  a  regular 
generalized  function  and 

h{y)=TG{y)-  (7) 

Following  [3]  and  [4],  we  define  now  the  product  of  expressions  Tp(y)  and 
TG(y)-  Let  H  be  an  (n  -|-  m)-dimensional  square  matrix  of  the  form 


H  :  = 
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where  M  is  a  matrix  of  dimension  nx  m,  all  elements  of  which  are  equal  to  zero, 

except  for  the  element  in  the  lower  left  corner,  equal  to  1,  and  Omxn  is  a  zero 

lil 

matrix  of  dimension  m  x  n.  Using  the  matrix  H,  we  define  quasi-derivatives 
(j  =  0, 1,  ...n  +  m)  of  a  given  function  y,  assuming 

and  quasidifferential  expression,  assuming 

rH(.y)  ■= 

It  is  well  known  that  the  domain  D{th)  of  the  expression  th  is  given  by 
D{th)  =  {y|y  G  D{tf)  and r^y  G  D{tg)}, 


and,  besides 

TH{y)  =  TG{TF{y)){-=  TGTF{y))  for  y  G  D{th){=  D{tgTf))- 

Now  we  consider  the  quasidifferential  equation 

rGijF{y))  =  Ay,  (8) 

where  A  is  a  complex  parameter.  We  use  the  symbol  y  to  denote  the  column  vector 
y  :=  {y'li  ,y^H  ^  transposition  symbol),  and  consider  the 

system  of  first  order  linear  differential  equations 

y'  =  (^^  +  A)y,  (9) 

where  A  :=  (Ajj)  is  a  square  matrix  of  dimension  (n  -|-  m),  whose  elements  are 
determined  by  the  equalities  A^+m,!  :=  A  and  Xij  :=  0  for  all  other  values  of  i 
and  j. 

The  conditions  (4)  and  (6)  imply  that  all  elements  of  the  matrix  H  are  locally 
integrable  on  functions.  Therefore,  the  conditions  (4)  and  (6)  ensure  the 
validity  of  the  existence  and  uniqueness  theorem  for  the  solution  of  the  Cauchy 
problem  for  system  (9)  stated  at  an  arbitrary  point  of  the  semi-axis  ii+.  On  the 
other  hand,  the  equation  (8)  is  equivalent  to  the  system  (9)  in  the  sense  that 
if  all  quasiderivatives  of  y  up  to  the  order  (n  -|-  m  —  1)  are  locally  absolutely 
continuous  on  and  y  satisfies  the  equation  (8)  almost  everywhere,  then  y 
is  a  solution  of  the  system  (9),  and  vice  versa,  if  y  is  a  solution  of  the  system 
(9),  then  its  coordinates  are  locally  absolutely  continuous  on  and  the  first 
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coordinate  y{=  satisfies  the  equation  (8)  almost  everywhere  on  Thus, 
the  conditions  (4)  and  (6)  imply  the  validity  of  the  existence  and  uniqueness 
theorem  for  the  solution  of  the  Cauchy  problem  for  the  equation  (8)  stated  at  an 
arbitrary  point  of  the  semi-axis  -R+. 

Formulas  (5)  and  (7)  allow  treating  the  differential  equation  (1)  with  distri¬ 
bution  coefficients  as  a  quasidifferential  equation  (8),  and  this,  in  turn,  allows 
asserting  that  the  conditions  (4)  and  (6)  imply  also  the  validity  of  the  existence 
and  uniqueness  theorem  for  the  solution  of  the  Cauchy  problem  for  the  equation 
(1)  stated  at  an  arbitrary  point  of  the  semi-axis 

3.  Main  results 

We  now  formulate  our  main  results.  The  symbol  o(l),  as  usual,  will  denote  a 
function  infinitely  small  as  x  — )•  -|-oo. 

Theorem  1.  Suppose  that,  the  number  A  in  the  equation  (1),  is  different  from 
the  product  anbm  CLnbm),  and  the  functions  pi,  p2,  ■  ■  ■  ,Pn  and  qi,  q2,  ■  ■  ■  ,qm 
satisfy  the  conditions 

n 

i=2 

and 

m 

+  (i  +  k2-gi|)X]l^il)  ^ 

i=2 

where  a  is  the  largest  multiplicity  of  a  root  of  the  polynomial 

d{z)  =  (z^  +  aiz^-^  +  ...  +  an)iz^  +  biz^-^  +  . . .  +  6^)  -  A. 

Then  the  equation  (1)  has  a  fundamental  system  of  solutions  yi,y2,  ■  ■  ■  ,yn+m, 
such  that  if  zi  is  a  root  of  the  polynomial  of  multiplicity  li  ^  a,  then  the 
equation  (1)  has  a  subsystem  {yj}  (j  =  1,  2, . . . ,  li)  of  fundamental  solutions  such 
that 

(x)  =  +  0(1)),  s  =  1,  2, . . . ,  n,  (10) 

y[n+^-l](^)  ^  ^  ^^^n+s-2  ^  ^  ^  ^ 

s  =  l,2,...,m, 

as  X  ^  +00. 

Another  subsystem  of  fundamental  solutions  {yff  {j  =  li  +  1, . .  .li  +  I2)  cor¬ 
responding  to  the  root  Z2  of  the  polynomial  ^{z)  of  multiplicity  I2  has  the  same 
asymptotics,  and  so  on. 
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Proof.  The  quasidifferential  equation  (8)  (and  so  the  differential  equation  (1)) 
is  equivalent  to  the  system  of  equations  (9).  We  write  this  system  in  the  form 

y' =  (A  +  i?(x))y.  (12) 

Here,  as  before,  y  is  the  column  vector  y  =  (y|J^ ,  ,  •  •  • ,  the  nonzero 

elements  Uij  of  (n  +  m)-dimensional  square  matrix  A  are  determined  by  the 
equalities 

Q'jj  +  l  —  j  —  f)2, ,71  +  771  1,  Onj  —  0,n+l—j  i  j  —  1)2,...,  77, 

(In+m,!  —  ®n+m,n+j  —  &m+l— j  —  1)2,  ...,777, 

and  the  nonzero  elements  Vij  of  the  matrix  function  R{x)  are  defined  by  the 
equalities 

Vn-lJ  =  -Pn+l-j,  j  =  1)  2,  ...77  -  1,  r^l  =  -{p2  -  Pi  -  ai)pn, 

rnj  =  Pn+2-j  -  {P2  -  Pi  -  ai)pn+l-j,  j  =  2,  3,  ...,  77-1,  Vnn  =  P2  -  Pi, 
'l"n+m—l,n+j  —  Qm+l—j,  j  —  1)  2,  ...777  1,  Tn.\-m,n+l  —  (^2  Ql  ^l)Pm) 

^n+m,n+j  —  Qm+2—j  (^2  ^l)ym+l— +  j  —  2, 3,. ..,777  1,  Tfi+m,n+m  —  Q2  ^1' 

Note  that  the  characteristic  polynomial  of  the  matrix  A  coincides  with  the 
polynomial  Let  zi  be  the  root  of  the  equation  ^{z)  =  0  of  multi¬ 

plicity  ^1,  i.e.  the  eigenvalue  of  the  matrix  A  of  algebraic  multiplicity  h.  The 
structure  of  the  matrix  A  is  such  that  ai^i+i  =  1,  ifl  <  i  <  77  +  777  —  1,  and 
ajj  =  0,  if  2  <  7  +  1  <  J  <77  +  777.  Therefore,  any  eigenvector  c  of  the  matrix  A 
corresponding  to  the  eigenvalue  zi  has  the  form 

c  =  {l,zi,z‘f,...,  z'f~^,Q{zi),  ziQ{zi), . . . ,  z^~^Q{zi))'^, 

where  Q{z)  =  +  aiz'^~^  +  . . .  +  a^. 

Thus,  only  one  eigenvector  corresponds  to  the  eigenvalue  zi  of  the  matrix  A, 
i.e.  the  geometric  multiplicity  of  zi  is  equal  to  1.  In  other  words,  each  eigenvalue 
of  the  matrix  A  is  associated  with  only  one  Jordan  block  in  its  canonical  form. 
Therefore,  the  dimension  of  the  Jordan  block  in  the  canonical  form  of  the  matrix 
A  of  greatest  dimension  coincides  with  the  multiplicity  of  the  eigenvalue  of  the 
matrix  A  of  greatest  multiplicity,  i.e.  with  the  number  a  from  the  condition  of 
Theorem  1.  Thus,  in  this  situation  the  system  of  equations  y'  =  Ay  has  solutions 
that  can  be  represented  as 

and  Vc  +  0(e"i  V”!),  j  =  1,  2, . . . , -  1. 
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In  addition,  it  follows  from  the  conditions  of  Theorem  1  that  all  elements  of  the 
matrix  belong  to  the  space  L^{R+).  Applying  the  statement  of  problem 

35  from  [5,  Ch.  Ill,  p.  120]  (see  also  [6,  Chapter  IV,  Remark  to  Theorem  4, 
p.  95]),  we  see  that  the  system  of  equations  (12)  has  a  subsystem  of  solutions  of 
{yj})  J  =  1)  2, . .  Ji,  representable  in  the  form 

Yj  =  +  o{l)),  j  =  1,  2, . . . , /i,  as  X  — > +00. 

A  further  look  into  the  relationship  between  the  solutions  of  the  equation 

(I)  and  the  solutions  of  the  system  (12),  shows  that  the  equation  (1)  has  a 
subsystem  of  solutions  {yj},  j  =  1,  2, ...  /i,  representable  in  the  forms  (10)  and 

(II)  as  X  — )■  +00. 

Repeating  these  arguments  for  an  eigenvalue  Z2  of  multiplicity  h,  we  find 
that  the  equation  (1)  has  a  subsystem  of  solutions  {yj},  j  =  h  +  1, ...  h  +  I2, 
with  an  asymptotic  form  (10)  and  (11)  (with  zi  replaced  by  Z2),  and  so  on. 
From  the  above  reasoning  it  also  follows  that,  having  considered  all  the  roots  of 
the  polynomial  5^(2:),  we  obtain  the  asymptotics  of  some  fundamental  system  of 
solutions  yi,y2, ... ,  yn+m  of  the  equation  (1).  The  theorem  is  proved.  ◄ 

Remark  1.  The  condition  A  /  of  Theorem  1  implies  that  the  number 

z  =  0  is  not  the  root  of  the  polynomial  ^{z).  This  allows  us  to  find  not  only  the 
main  term  of  the  asymptotics  of  the  solutions  of  the  equation  (1),  but  also  their 
quasiderivatives  up  to  the  {n  +  m  —  l)th  order.  In  particular,  asymptotic  formulas 
for  solutions  of  the  equation  (1)  (formulas  (10)  with  s  =  1)  can  be  differentiated 
up  to  the  order  (n  —  2).  If  X  =  anbm  and  the  number  z  =  is  the  root  of  the 
polynomial  S^(z)  of  multiplicity  I,  then  the  method  of  the  proof  of  Theorem  1  allow 
us  to  state  that  the  equation  (1)  has  solutions  of  the  form 

yj  =x^~^{l  +  o{l)),  j  =  l,2,...,l, 

but  for  quasiderivatives  of  these  solutions  one  can  only  say  that 

y^j\x)  =  x^~^o{l),  s  =  1,2, . . .  ,n  —  l,n  +  1, . . .  ,m  +  n  —  1, 

yf\x)  =x^"^(a„  +  o(l)). 

Using  more  subtle  methods  for  first-order  linear  differential  systems,  one  can 
also  study  the  situation  when  A  =  In  particular,  the  following  theorem 

holds. 

Theorem  2.  Suppose  that,  in  the  equation  (1), 

ai  =  02  =  ...  =  On  =  bi  =  b2  =  ...  =  bm  =  A  =  0 
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and  the  functions  pi,p2,  ■  ■  ■  ,Pn  o,nd  qi,  q2,  ■  ■  ■ ,  Qm  satisfy  the  conditions 

n 

\pi\  +  {1  +  x\p2  -  Pi\)^x^~^\pj\  G  L\R+) 

J=2 

and 

m 

\qi\  +  (1  +  x\q2  -  qi\)^x^~^\qj\  G  L^{R+). 

i=2 

Then  the  equation  (1)  has  a  fundamental  system  of  solutions  {yj}  (j  =  1,  2, . . . ,  n+ 
m)  such  that 


(x) 


(j  -  1  -  s)\ 
o(l), 


(1  +  o(l)), 


if  s  =  0,1,..., j  -  1, 

if  s  =  j,j  +  l,...,n  +  m-l. 


(13) 


as  X  ^  +00. 

Proof.  As  we  already  noted  in  the  proof  of  Theorem  1,  the  equation  (8)  (and 
so  (1))  is  equivalent  to  the  system  (12).  Moreover,  it  follows  from  the  conditions 
of  Theorem  2  that  the  matrix  A  has  a  Jordan  form  containing  one  Jordan  block 
with  zero  eigenvalue  of  multiplicity  (n  +  m).  By  the  symbol  D  we  denote  a 
diagonal  matrix 

D  :=  diag{l,  x,. . . ,  x”, . . . , 

The  calculations  show  that  the  nonzero  elements  fij  of  the  matrix  R  := 
DRD~^  are  determined  by  the  equalities 

rn-i,j  =  -x'^~^~^Pn+l-j,  j  =  1,2,. . .  ,n  -  I, 

rn,j  =  x"‘~^{pn+2-j  +  {pi  -  P2)Pn+l-j),  j  =  2,  3,  .  .  .  ,  U  -  1, 
rn,l  =  X"'~^{pi  -p2)Pn,  fn,n  =  P2  -  Pi, 
i’n+m— l,n+j  —  X  ^  qm+l—j,  j  —  1)2,  ...771  1, 

i'n+m,n+j  —  X  ^  (^qm+2—j  T  {Q1  Q2')Qm+l—j) ,  j  —  2,  3,  ...,  171  1, 

^n+m,n+l  ~  X  [qi  q2)Pm,  Xn+m,n+m  —  Q2  91- 

From  these  formulas  and  from  the  conditions  of  Theorem  2  it  follows  that  all 
elements  of  the  matrix  R  belong  to  the  space  L^(i?_|_).  Thus,  the  coefficients  of 
the  system  (12)  satisfy  all  the  conditions  of  Theorem  1.10.1  from  [7].  Applying 
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the  above  theorem  again,  we  see  that  the  system  (12)  has  solutions  {yj},  j  = 
1,  2, n  +  m,  of  the  form 


yj  = 


U  -  1)! 


A  further  look  into  the  relationship  between  the  solutions  of  the  equation  (1) 
and  the  solutions  of  the  system  (12)  shows  that  the  equation  (1)  has  a  system  of 
solutions  Hj,  j  =  1,2, . .  .n  +  m,  representable  in  the  form  (13)  as  x  — )•  +oo.  The 
theorem  is  proved.  ◄ 

It  will  be  interesting  to  join  the  method  of  this  paper  with  method  of  paper 
[8]  for  investigation  of  more  complicated  equations. 
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